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ABSTRACT 


A system for measuring unsteady pressures in Tlow fields 
employing a remote transducer and thin plastic pressure transmitting 
lines has been designed and built. 

The static and dynamic characteristics of the system were 
determined experimentally, and the results were found to agree well 
with a theoretical model. 

The measuring system was subsequently integrated into a two- 
dimensional wind tunnel mode! consisting of a symmetrical airfoil with 
a plain oscillating flap. The flap was harmonically oscillated, and 
the measuring system was used to determine both the steady and unsteady 
pressures at a point near the flap hinge line. The static and dynamic 
results were then compared to those obtained using thin airfoil theory 


and found to agree well when corrected for thickness effects. 
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TABLE OF SYMBOLS 


Latin Symbols 


ay Mean velocity of sound 

e Airfoil chord 

Cy Specific heat at constant pressure 

AC, Pressure coefficient across airfoil (non-dimensional) = = 
C Specific heat at constant volume 

D Tube inside diameter 

g Gravity constant 

1 Imaginary unit = [-1 

J, Bessel function of first kind of order n 
k Polytropic constant for volumes 

k Non-dimensional frequency = oT 

te Length of flexible tube 

p Static’ pressure, including steady and unsteady terms 
PO Mean pressure 

P Amplitude of pressure disturbance 

P, Driving pressure 

Py Measured pressure (transducer sensed) 

Pr _ Pressure on airfoil lower surface 

Eu Pressure on airfoil upper surface 

2s Prandt] number = Mave 

q Dynamic pressure = y pu? 

r Radial coordinate 

R Tube radius 

R Gas constant 

ip Time 


—t 


To 


(a | 


Combined temperature inside tube, steady plus unsteady component 
Mean temperature inside tube 

Amplitude of temperature disturbance inside tube 

Complex velocity component in axial direction inside tube 
Amplitude of velocity component in axial direction inside tube 
Velocity of fluid meeting the airfoil 

Complex velocity component in radial direction inside tube 
Amplitude of velocity component in radial direction inside tube 
Volume of pressure transducer 


Tube volume 


Greek Symbols 


x 


B 


nf 


C as 


Di 


Wave number (Equation 6) 

Location of chordal station in transformed coordinates 
Ratio of specific heats -2 

Control surface deflection, positive deflection trailing edge down 
Phase angle 

Coefficient of thermal conductivity 

Absolute fluid viscosity 

Location of hinge line station in transformed coordinates 
Frequency 

Fluid density 

Mean fluid density 

Amplitude of fluid density disturbance 


Time constant 


CHAPTER I 


INTRODUCTION 


There is a clear need for a practical and inexpensive tech- 
nique to measure unsteady pressures on experimental models. The 
character of the pressures being measured may be either well defined 
and harmonic, as might be found on an oscillating airfoil, or random 
natured as found in the separated flow region of flow about a bluff 
body. Although progress has been made in the fabrication of small 
light weight transducers which may be installed directly in a model, 
the problem remains of accounting for the model dynamic effects on 
the transducer response. Two other factors limit the use of this 
type of system: space limitations, such as those encountered when 
working with a thin airfoil, and the important factor of cost. 

It is a commonly accepted technique to use long thin pressure 
leads’ to connect a static orifice to a centrally located pressure 
sensing instrument such as a manometer board, a micromanometer, or 
an electrical pressure transducer for the purpose of measuring average 
pressure. Although it is apparently not well known in this country, 
a suitably sized pressure lead used in conjunction with a properly 
selected transducer cavity volume will also provide adequate dynamic 
response traits so as to be useful for measuring unsteady pressure. 

Bergh [1] at the National Aeronautical and Astronautical 
Research Institute (N.L.R.) in Amsterdam, Holland, has demonstrated 
a technique that employs long thin pressure leads for the measure- 
ment of unsteady pressures. He has presented measurements of the 
in-phase and out-of-phase pressure distributions on harmonically 


oscillating Tee-tail configurations in’ a sub-sonic wind tunnel. The 


li 


technique employs multiple pressure leads, having identical dynamic 
traits, consecutively sampled in a scanning valve arrangement by a 
Single pressure transducer system. Similar methods’ have been employed 
by Laschka [2] on a moving rigid model at the Entwicklungsrings Sud in 
Munich,’ Germany, with good agreement between theory and experiment. 
Although employed in European laboratories, this’ technique has appar- 
ently not been used to any extent in this country. 

It was the purpose of the investigation described herein to 
confirm the results of Bergh [1,3] by direct frequency response 
measurements upon a home built pressure measuring system and then 
apply the technique of measurement to a typical unsteady flow environ- 
ment. For the latter purpose, a two-dimensional airfoil model with 
an oscillating plain flap was designed and instrumented with pressure 
orifices just ahead of the hinge line. The results obtained experi- 
mentally are then compared to those obtained using thin airfoil theory. 

Because of a concern for using the technique in environments 
other than harmonic flow, the indicial response traits of the system 


have also been determined. 
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CHAPTER I] 


EXPERIMENTAL EQUIPMENT AND PROCEDURES 


IT. EQUIPMENT 
General 

The work described here was conducted in the Dynamics 
Laboratory and the 3.5 foot by 5.0 foot wind tunnel of the Naval 
Postgraduate School's Aeronautical Engineering Laboratories. 

The dynamics laboratory provided the work area for construc- 
tion and calibration of the locally prepared pressure pickups. The 
wind tunnel provided a place to demonstrate the pressure measuring 
method on a practical model, in this case, a two-dimensional flapped 


airfoil. 


Pressure Transducers 

The pressure pickups employed in the study described here 
were fabricated from a local design shown in Figure 1. The transducer 
assembly consisted of a Bentley Nevada Corporation Model D-252 Distance 
Detector system, a brass sensing head holder and a 0.003 inch thick 
annealed brass’ diaphragm. The Distance Detector system was composed 
of a D-252 Distance Detector, a 15 volt direct current regulated power 
supply and a H-1-022-4 sensing head. 

The distance detector system is a reluctance gage which senses 
changes in reluctance produced by changes in air gap between the sensing 
head coil face and the surface of a conducting material. Brass was 
chosen as the diaphragm material in this case due to its availability 
and ease of fabrication. Other materials which could equally wel] 
have been used as a reference include gold, silver, platinum, copper 


and aluminum. 
is 
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The annealed brass diaphragm was press formed, cut to size 
and sweat soldered to the face of the brass sensing head holder. 
Brass was selected for the sensing head holder in order to reduce 
differential expansion effects caused by thermal changes. After the 
diaphragm was attached to the holder, the sensing head was inserted 
in the holder as may be seen in Figures 1 and 2 which show the pres- 
sure transducer. 

During the initial stages of transducer construction a dia- 
phragm thickness of 0.003 inches was chosen. The choice of diaphragm 
thickness was based upon the need for a high sensitivity in the pres- 
Sure range, referred to ambient, of plus or minus 1.0 pound per square 
inch (p.s.i.). The transducer constructed for use here is similar to 


one used by Schmidt [4] in an investigation of fluctuating air loads. 


Calibration Chamber and Transducer Holders 

The calibration chamber was designed to determine the static 
and dynamic characteristics of the pressure transducer systems. It 
was constructed of an aluminum block 4 inches long with a 1.5 inch 
diameter cylindrical chamber in the center and is shown in Figure 3. 
The chamber had both ends faced and counter bored to accept either a 
transducer or an adapter plate used to connect the chamber to a pres- 
sure line.” Two additional openings were provided in the calibration 
chamber: one to accept an acoustical driver unit, in this case a 
University Sound ID-75 Driver Unit, and the other a pressure port to 
enable the mean static pressure in the chamber to be varied. The 
acoustic driver unit was. placed. at the midpoint of the acoustic cavity 
length so that unsteady pressure pulses could be introduced symmetri- 


cally into the cavity. 





FIGURE 2 PRESSURE TRANSDUCER 
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During static calibration both transducers were installed 
directly in the chamber, while for the dynamic ca’*bration the adapter 
plate was placed over one end of the chamber and tne remote transducer 
placed in a suitable hoider. The noider was then connected to the 


calibration chamber using a flexibie Tygon tube 


Instrumentation 

The output from a pressure transducer was ted into an Astro- 
Data Model D-885 amplifier. The output from the Astro-Data amplifier 
was read on a United Systems Corporation. Digitec, digitai voitmeter 
for the static calibration and on a Ballentine Laboratories Mode! 320A 
True root mean square (R.M.S) meter during dynamic calibration. During 
dynamic calibration the signal from a Hewiett Packard Model 200CD wide 
range oscillator was amplified by a Genera’ Radio Company Type 1233-A 
power amplifier and used to actuate the acoustical driver unit, and 
an AD-YU Type 405 precision phase meter was used to measure the phase 
difference between the output voltages of the two transducers. The 


instrumentation is shown schematicaily in rigures 4 and 5. 


Wind Tunnel Model 

The wind tunnel model used was a modified NACA 64 -Ole airtoil 
with a 20 per cent chord, radius nosed, plain flap. The model had a 
span of 22 inches with a constant chord of i2 ‘nches The tips of the 
model were terminated in end plates in order to simulate two-dimensional 
flow traits. Modification of the airfoil section was achieved by 
Straight line fairing the surface from the 60 per cent*chord position 
to the trailing edge to remove the cusp. Two flush mounted static 
pressure taps were installed on the upper and lower surfaces of the 


airfoil 0.6 inches ahead of the flap hinge line (0.75 per cent chord 


18 


Pressure 


Sogrce ee —— 
| 


Direct Reading 
Micromanometer 


Pressure Transducer 
alibration Pl 
Chamber 


Bentley Model D252 
Distance Detectors 













-15 Volt D.C. Regulated Power Supply 


Bias Voltage 


Astro-Data Model D885 
t Amplifiers f 


Digitec 
D.C. Voltmeter 


FIGURE 4 STATIC CALIBRATION INSTRUMENTATION 


19 


Hewlett Packard Model 200CD 
Wide Range Oscillator 


General Radio Company Type 1233A 
Power Amplifier 


~ = University Sound 
ID75 Acoustical Driver 


ingen Pressure Transducer 


aq Calibration! 
Chamber |_| 











Bentley Model D252 
Distance Detectors 


-15 Volt D.C. Regulated Power Supply 


{1 E- 


mn Bias Voltage Gis 


L 


—— ew 
<= 


! 


: Astro-Data Model D885 
Y Amplifiers 





Dual Beam AD-YU Type 405 Ballentine Model 320A 
Oscilliscope Precision Phase True R.M.S. Meter 
Meter 


FIGURE 5 DYNAMIC CALIBRATION INSTRUMENTATION 


20 


position). These taps were connected to the pressure transducers with 
flexible Tygon tubing. The tubing used was a soft plastic with a 
0.0465 inch inside diameter and a 0.125 inch outside diameter and had 
been previously calibrated relative to its unsteady pressure frequency 
response traits. The airfoil section is shown in Figure 6. Figures 7 
through 10 are photographs of the wind-tunnel model employing the pres- 
Sure measuring system. The plain flap was of light weight balsa 
construction glued to a 0.125 inch diameter steel drill rod. The 
deflection angle was sensed by a potentiometer located at the left- 
hand end of the model. A small steel] tab installed in the center of 
the trailing edge of the flap and connected to a constant speed motor 


with steel piano wire provided actuation of the flap for dynamic tests. 


II. CALIBRATION PROCEDURES 

Static’ Calibration 

Static response of the two pressure transducers used in this 
study was determined by installing the transducers in the calibration 
chamber, applying a known pressure, measured with a direct reading 
micromanometer, to the chamber and’ measuring the output voltage. The 
air gap between the sensing head and the brass diaphragm in the trans- 
ducer was adjusted during this procedure to obtain a transducer output 
of 0.28 volts per p.s.i. which was linear in the’ pressure range of plus 
or minus 1.2 pounds per square inch gage (p.s:.i.g.). A pressure change 
of 1.0 p.s.i. deflected the diaphragm approximately 0.00065 inches 
which corresponds to a deflection of less than one third the diaphragm . 
thickness for the desired full pressure range. It was felt that an 
upper limit of one third the diaphragm thickness was the maximum deflec- 


tion allowable in order to maintain a linear pressure calibration. 
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To obtain the dynamic characteristics of the pressure trans- 
ducers and the proposed measuring system, the installation shown in 
Figure 5 was used. The initial runs were conducted with both trans- 
ducers mounted in the calibration chamber. This configuration allowed 
a comparison to be made of the characteristics of the two transducers. 
The symmetry of the acoustical wave propagation was also verified by 
this test. 

In order to obtain the dynamic response of the system to a 
remote pressure pickup, one transducer was connected to the calibra- 
tion chamber by a flexible plastic tube as described above. Two dif- 
ferent diameters of tubing were used, and the cavity volume which is 
formed in front of the transducer diaphragm face was varied, using 
ringlike spacers of varying thickness. 

For dynamic calibration runs a fluctuating chamber pressure 
was produced by the acoustical driver unit. The values of phase 
difference and the R.M.S. value of pressure for both transducers were 
recorded at frequency increments of 5 Hertz from 10 to 400 Hertz. 

Using the data recorded in each run a system dynamic gain, 
defined as Po/P. > was calculated. The values of calculated gain and 
measured phase angle were then plotted as a function of frequency to 


give a calibration curve for the system. 
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CHAPTER III 
RESULTS AND DISCUSSION 


Static Calibration 

Figure 11 is the measured static calibration curve for the 
pressure transducer. The transducer output was amplified by a factor 
of ten, and the voltage output as a function of pressure was plotted 
to give a calibration curve. The curve was linear in the range plus 
or minus 1.0 p.s.i.g. This curve was checked repeatedly during a 


period of three months, and the variation was less than one per cent. 


Dynamic Calibration 


Dynamic calibration began by checking the dynamic response of 
the two transducers when subjected to the same fluctuating pressure. 
It was found that both systems had identical responses trom 0 to 2500 
hertz with the exception of some erratic behavior with respect to pres- 
sure signal wave form between 1500 and 2000 hertz, but this frequency 
range is near the resonant frequency of the calibration chamber cavity 


which, for standard atmospheric conditions, iS given by: 


a 
e s) 
gmsescnance = (twice the cavity length) ° (1) 


and equals 1835 hertz. 
The harmonically varying pressure at the input end of a tube 
is given by: 
“ae iwt 
P= Ee +P, e ; (2) 
and at the tube exit into the transducer cavity volume is: 
P,= Po +P, el (werd) (3) 
Therefore, the dynamic characteristics of the pressure measuring system 


P 

are defined by the amplitude ratio, py and the relative phase shift, 6. 
I 
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In order to determine the effects of very short rigid pres- 
Sure. lines, an adapter plate containing a 0.020 inch diameter hole was 
used to separate one transducer from the chamber . The adapter had shim 
inserts which allowed the distance that the transducer cavity was sepa- 
rated from the calibration chamber to be varied from 0.188 inches to 
0.288 inches. It was found that the oscillatory signal was trans- 
mitted to the remote transducer almost unattenuated from 0 to ap- 
proximately 100 hertz with a phase lag of 6.0 degrees or less, with 
the measured signal lagging the applied signal. This result confirms 
results obtained by Schmidt [4] and stated by Berg and Tijdeman [3]. 


P 
Figure 12 shows a typical plot of p~ versus w. This case is for an 


4 





orifice diameter of 0.020 inches, ltength of 0.288 inches and trans- 
ducer cavity volume of 1.454 x 10° cubic inches. The transducer 
cavity was formed using a Teflon washer 0.038 inches thick with a 
0.60 inch diameter hole in the center. Other volumes were obtained 
by varying the thickness of the washer. 

The next phase of determining the behavior of oscillatory 
pressure transmission through a thin tube consisted of moving the 
second transducer away from the calibration chamber to an external 
holder connected to the calibration chamber by a plastic tube as 
described above. Three parameters were studied during this phase: 
tube diameter, tube length and transducer volume. Figures 13 and 
14 show the effects of varying the individual parameters on the 
system response. It is seen that increasing the length of tube or 
the transducer cavity volume decreases the ampiitude ratio and shifts 
the peaks to a lower frequency, while increasing tube diameter in- 


creases amplitude ratio and moves the peaks to the higher frequency. 
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PRESSURE RATIO, Pa/P, 


SYSTEM CHARACTERISTICS 






| | A - L = 36.45 in., V=1.454 x 102 in., 
a D = 0.065 in. : 
—T— f= L = 36.45 in., Vj=2.956 x 10° in., 
: D = 0.065 in. fee 

| @- = 40.0 in., V=1.454 x 10 in., 

DD = 0.0465 in. 7 
@- L = 93.45 in., V=1.454 TG. in. 

D = 0.065 in. 


50 T00 150 200 250 
FREQUENCY, HERTZ : 


FIGURE 13 EFFECT OF SYSTEM GEOMETRY ON PRESSURE 
RESPONSE 


32 


| 
: 
: 


PHASE SHIFT, DEGREES 


SYSTEM CHARACTERISTICS 





#2 3 
L = 36.45 in., V.=1.454 x 10 in., 
D = 0.065 in. 2 3 
O- L = 36.45 in., V=2.956 x 10 in., 
-D = 0.065 in., — 
480 | - ®- Ll = 40 in., V=1.454 x 10 in, 
| D = 0.0465 in. ane 
®- L = 93.45 in., V=1.454 x 10 in., 
| D = 0.065 in. 
400 | —— — 
| 
ye 
320 
yp | 
240 a 
a | cr : us A = 
eA x Pies ee 
10} - / —— geERr 1- 
Mie 
iG os 
/ ot ar | 
g a Bageg fe 
comments | 
80 : “et ’ Jl ‘ ee —— : _- 
rT" 
oe 
ae We 
vad 4%" 
0 50 100 150 200 250 


FREQUENCY, HERTZ 


FIGURE 14 EFFECT OF SYSTEM GEOMETRY ON PHASE SHIFT 
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The relative phase shift was only slightly changed by varying cavity 
volume and tube radius, but a much greater phase shift was measured 
for a given frequency when the tube length was increased. 

The removal of the transducer trom the calibration chamber and 
its installation in an external holder now simulates the working system. 
The unsteady pressure signal is now transmitted from the source (the 
calibration chamber) to the measuring device (the transducer) through 
the thin plastic tube. 

In choosing a working system it was feit that a smooth transfer 
function was necessary for practical use Therefore, if a combination 
of parameters produced sharp resonance peaks during dynamic calibra- 
tion, the system, although we!! defined, was considered an unusable 
configuration. 

One approach that was used in an attempt to smooth out reso- 
nance was the placing of a tee in the line between the pressure source 
and the transducer. The third leg of the tee was fitted with a tube 
which was to act as an acoustic damper. This idea proved unusable as 
sharper resonance peaks were developed no matter what the length of 
the damping tube, which was varied from 10 per cent to 100 per cent of 
the length of the pressure lead. An even more disturbing effect from 
the system with a tee was the erratic phase shift produced, possibly 
due to the different values of resonance associated with the different 


tube’ lengths. 


Wind Tunnel System Calibration 


The pressure measuring system chosen for use in. the wind tunnel 
was composed of 40 inches of 0.0465 inch inside diameter plastic tubing 


connected to a transducer volume of 2.965 x 10° cubic inches. The cali- 


bration curve for this system is shown in Figure 15. 
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PHASE SHIFT, DEGREES 


Repeatability of the calibration run was found to be good, as 
the variation during runs was less than one per tent. Table I sum- 
marizes the experimenta! calibration vaiues of phase angle and ampli- 
tude ratio for each frequency that generated the calibration curve 
applicable to the configuration used in the wind-tunnei measuring 
sys tem. 

A theoretical model to predict the response of a remote pres- 
Sure measuring system to an oscillating pressure has been proposed by 
Bergh [1,2]. In Bergh‘s model (See Appendix A) the amplitude ratio 


is given by: 


P V J,{4) 
7 wt YW LE Yola) | = 
me [cosh(oL) + KV, 9 ortay aeel (4) 
| o 
where 
1 
hela. * 25 -1 Je (oP,2) 3 
Ww 2 4 ‘~] pr“) 4% 
Stele i 2 
W 1 
a = iFR(PO*)* (6) 


In order to compare the theoretica! mode! with the curves 
obtained for the actual system, two theoretical points were calcu- 
lated. The calculations’ were laborious, and a digitai computer 
program would have been desirable for calculating more points and 
establishing a theoretical curve. The calculated points are shown 
on the calibration curve for the system (symbols ©, 0). 

The calculated points are in reasonabie agreement with the 
calibration curve. The discrepancies which are present may be at- 
tributed mainly to the uncertainties in the measured physical prop- 


erties of the system. The smail diameters and thickness which must 
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TABLE I 
DYNAMIC CALIBRATION DATA FOR PRESSURE MEASURING SYSTEM WITH 


L = 40 inches Ne = 2.956 x liom Cubic inches 
D = 0.0465 inches 


Frequency Pressure Ratio Phase Shift 

(Hertz) (P,/P. ) (degrees) 
oo P, lags P 

ome 
10 1.06 14.5 
15 1.09 aa 5 
20 kes S| .0 
25 edligs 40.0 
30 1.219 49.5 
35 1.244 61.5 
40 1.244 ciple 
45 lee08 86.5 
50 1.144 99.0 
55 ].063 109.5 
60 0.969 102.0 
65 .888 128.0 
70 813 | S90 
75 750 143.0 
80 . 700 149.5 
85 .650 ba0 
90 Ole 162.0 
95 .588 167.0 
100 .566 17a0 
100 560 1¢7.0 
110 .540 189.0 
120 BS) Sis 198.0 
130 .540 209.0 
140 .560 aC) 
150 565 238.0 
160 sos 2 Sault 
170 ns 6 5 ra) Sie 
180 ~ 485 291.0 
190 .430 306.0 
200 ged 319.0 
ZN - 348 oi) a0 
220 .318 341.0 
230 295 301.0 
240 285 364.0 
250 280 50/0 
260 .278 378.0 
270 .278 389.0 
280 .280 403.0 
290 0.280 417.0 
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be measured could easily lead to inaccuracies in calculated data. For 
example, the inside diameter of the plastic tubing could only be meas- 
ured using a wire drill as a plug gage Ths at best gave a rough mean 
inside diameter. The thickness of the Tefion used to produce the 
transducer volumes was found to vary in thickness by three to four 
thousandths of an inch. Thus the thickness was a mean value. 

Bergh and Tijdeman [3] present a number ot caiibsation curves 
for various pressure measuring systems Although no attempt was made 
to exactly duplicate any one of the systeins investigated by Bergh and 
Tijdeman, the same que‘itative System responses were noted in this 
Study. In the present study the purpose was to prove experimenta! ly 
that a workable pressure sénsing System couid be constructed and cali- 
brated with a high degree of reitability That the system was reliable 
is demonstrated by the fact that during the experiments both the static 
and dynamic calibration curves were frequentiy checked and found to 
vary less than one per cent. 

In order to extend the system capability to measure random 
pressure fluctuations the indicial response of the system was measured. 
The calibration chamber pressure was abruptly changed from one p.s.i.g. 
to atmospheric. and the system response: was displayed on an oscilliscope. 
It was found that the system time constant, 7, where: 

+ ='At (for pressure drop from P, to atmospheric), (7) 
was thirteen milliseconds with the transducer instailed directly in 
the calibration chamber and fourteen miiliseconds with the transducer 
connected to the chamber by a 40 0 inch tube, demonstrating that the 
effect’ of the connecting tube on system response is indeed smali. A 


“dime store" balloon connected to the calibration chamber through a 
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0.25 inch diameter tube was exploded in order to produce the abrupt 


pressure Change. 


An Application of the Pressure Measuring System 


Once a dynamic calibration curve had been obtained, the meas- 
uring system was integrated into the wind-tunnel model. 

A zero frequency measurement was performed to determine the 
Slope of the static pressure coefficient between the top and bottom 
surface per degree of contro! surface deflection. The pressure coef- 


ficient 1s given by: 
Pi - Py 


(8) 


ACp = 


The position of the pressure taps on the airfoil were at the 75 per 
cent chord position which was 5 per cent chord ahead of the control 
surface hinge line. To verify that the static pressure measurement 
was correct, pressure taps similar to those connected to the pressure 
transducers were installed at the same chord location and connected 
to’ a manometer board. The pressure change across the airfoil measured 
in this manner was found to be the same as that measured by the pres- 
Sure transducers. The value of the slope of the static pressure coef- 


ficient curve per degree of flap deflection was found to be 


dé static 0.0ay peradegnec: 


Thin airfoil theory, Appendix B, indicates a value for this 
coefficient of 0 059 per degree, which is larger than that actually 
obtained. Etkin [5], however, discioses that thickness ratio has an 
effect upon the hinge moments of a flapped airfoil, and it would be 
expected that the static pressures near the hinge line of the control 


Surface would be similarly affected The hinge moment value given in 


oe 


[5] for'a 12 per cent thick airfoil is 83.5 per cent of the zero 
thickness value. If this correction were applied to the calculated 
Slope of the pressure coefficient, the value would become 0 049 per 
degree which agrees well with the measured value. 

In order to determine the effect of control surface gap on 
the static pressure readings, a zero frequency. test was conducted 
for varying surface deflections with the gap completely sealed. The 
Static pressures obtained during this test were identical to those 
obtained with the gap unsealed, and it can be assumed there was little 
or no flow leakage through the control surface gap. 

Due to mechanical limitations it was possible to run at only 
a single frequency to check the dynamic characteristics of the pres- 
sure measuring system. The control surface was oscillated at 28.5 
hertz using a constant speed electric motor. This limitation of a 
single frequency would be a deficiency in an experiment if the goal 
were other than the operational proof of the system. The R.M.S. values 
of control surface deflection and upper and lower surface pressures 
were recorded for eight runs with a tunnel dynamic pressure of 24.18 
pounds per square foot. The measured pressures were corrected for 
transducer’ system dynamic gain according to Figure 15,°'the amp? itude 
ratio, Bh for 28.5 hertz being equal to 1.20. Using the data obtained, : 
a dynamic value for the pressure difference across the wing was deter- . 
mined. The average value of the slope of the pressure coefficient per 


degree was 


dep dynamic 7 0.045 per degree. 


Figure 16 is a simultaneous recording made of the output of both pres- 


Sure transducers and the control surface position on a Brush Instruments 
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The following R.M.S. values of pressure and control 
surface deflection were recorded with tunnel] dynamic 
pressure, q, equal to 24.18 pounds per square foot 
and flap frequency, a, equal to 28.5 Hertz. 


Lower Surface Pressure = 0.020 p.s.1i.g. 
Upper Surface Pressure = 0.020 p.s.i.g. 
Control Surface Deflection = 4.99 degrees 


FIGURE 16 OSCILLOGRAPH RECORDING OF FLAP DEFLECTION 
AND PRESSURE OSCILLATION ON THE AIRFOIL SURFACE 
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Ultralinear Oscillograph. From left to right the traces are of pres- 
sure fluctuation on the lower surface, the upper surface. and control 
surface displacement. The traces are annotated to indicate maximum 
control surface travel, which was approximately 7.0 degrees in either 
direction, and maximum positive and negative pressures. The oscil- 
lograph traces were used as a check on the measured oscillation fre- 
quency and to monitor the behavior of the transducers and control 
Surface driver unit. The values of pressure used in obtaining pres- 
sure coefficients were those obtained from 2.M.S. readings. Kussner 
and Schwarz [6] present a method for determining the pressure coef- 
ficient across an airfoil with an oscillating plain flap (See Appendix 
B). The pressure coefficient is shown to be a function of the non- 
dimensional frequency, k, of the airfoil and is given by: 


ACp(8,t) = - 6 e | trRe(k,8,W) + iIm(k,8,) J (9) 
where 

cos B= - x/c. (10) 
This equation comes from the paper by Kussner and Schwarz [6], may 
also be found in Fung [7], and is more completely defined in Appendix 
B. The results of this equation are shown on Figure 17, where the 
pressure attenuation and phase angle are plotted versus non-dimensional 
frequency, k, for a value of 8 corresponding to the 75 per cent chord 
position when the plain flap hinge iine is at the 80 per cent chord 
position. The value of non-dimensional frequency which corresponds to 
the measured wind tunnel data is 0.77, which gives a pressure attenua- 
tion of 0.93. Using the value of the slope of the static pressure 


coefficient curve calculated from thin airfoil theory and corrected 


for finite thickness and applying the attenuation factor of 0.93, the 
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FIGURE 17 PRESSURE ATTENUATION AND PHASE 
ANGLE VERSUS NON-DIMENSIONAL FREQUENCY 
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calculated value of the slope of the dynamic curve becomes 0.046 per 
degree which agrees well with the experimental results. The pressure 
peak at the 75 per cent chord position should lead the peak flap oscil- 
lation theoretically by 27.5 degrees. No attempt, however, was made to 
measure this phase difference. 

A final check on the system calibration was performed when the 
work in the wind tunnel was complete. The calibration curve, Figure 15, 


was found to be confirmed to within one per cent. 
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CHAPTER IV 
CONCLUSIONS AND RECOMMENDATIONS 


I. CONCLUSIONS 

Calibration showed the static and dynamic properties of the 
pressure transducer system very well behaved. Both steady and unsteady 
static pressures were easily measured, and the measurements were 
repeatable to a high degree of accuracy. The correlation between 

theory and experiment could have been better, but precise measure- 
Pci of the physical parameters were difficult to obtain. Indeed, 
for a large number of tubes in an installation employing only one 
transducer and a scanner valve type arrangement, it would be advan- 
tageous to install a second transducer internally to the model. This 
second transducer would allow more flexibility in the use of the 
measuring system. A continuous calibration of the system would be 
available by simultaneously measuring the pressure at a given station 
with both transducers. 

The system is very easily installed in a model. Weight and 
cost wise the advantages are great. Unsteady pressure loads on an 
oscillating model are easily obtained without the necessity of cor- 
recting readings for the acceleration response of the individually 
mounted internal transducers. For these reasons, the system would be 
good for determining dynamic loadings for flutter analysis without the 
need to build an elaborate dynamic model with many small and costly 
model installed transducers. 

The experimental results obtained during the wind tunnel tests 


compare well with those obtained using thin airfoil theory. 
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Use of this method of unsteady pressure measurements is not 
limited to cases of purely harmonic flow. The calibration curve 
generated for the measuring system is in essence a system transfer 
function. With the transfer function known, the measuring technique 
may be used to measure stochastic pressure variations such as those 
encountered in regions of separated flow. For this type application, 
instead of measuring the transducer output directly, a recording of 
the transducer output could be made on magnetic tape and a power 
Spectral analysis performed on the recorded signal to. determine the 


character of the random pressure variation. 


TI. RECOMMENDATIONS 

Pressure measuring systems of this type have not been used 
before at the Naval Postgraduate Schoo!. Neither have they been used, 
as far as it is known, to any great extent in this country. The 
system constructed for this study has proven practical but could be 
improved upon in a number of ways. 

In the area of transducer construction a smalier sensing head 
could be utilized to construct a smaller transducer. The instrumenta- 
tion setup seems’ to be sound, but some transient high frequency noise, 
although only in the 2.0 to 3.0 millivolt range, was a bother and could 
be eliminated with some work on the circuit wiring. 

Bergh's [1] theoretical analysis appears accurate. If a 
digital computer program were written to use this analysis in a parame- 
tric study of the physical variables, it might be much faster to choose 
a suitable combination of these variables to fit a particular applica- 


tion. 
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This measuring system should be integrated into a larger and 
more complex model in conjunction with a multi-channeled scanning 


valve arrangement, to explore the full capability it seems to have. 
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APPENDIX A 


THEORY FOR DYNAMIC RESPONSE 
OF A THIN TUBE TO OSCILLATING PRESSURE 


The treatment presented here is due to Bergh [1] and is 
concerned with the solution for the pressure response into a trans- 
ducer cavity when a long slender tube is used to connect the trans- 
ducer cavity to the unsteady pressure source. No correction is made 
for tube wall elasticity as might be envisioned in a biological 
environment since it was experimentally confirmed that the difference 
in dynamic response between using a metal tube or a plastic tube was 
negligible for the configurations considered. 

The unsteady terms are postulated as harmonic perturbations 
Superimposed upon the mean average values in the following expres- 


sions for the dominant variables: 


Pressure: a 2 elms 
Density: al aan a 
TK iwt 
Temperature: 1 ie sa Ss (A.01) 
Axial Velocity: u= uy + u g wt 
Radial Velocity: v = ie + Vv sien 


The above terms are considered in the polar coordinate form 
due to the inherent symmetry about the tube axis. Typically, the 
pressure perturbation amplitude, p, in equation (A.01) would be func- 
tionally expressed as p(x,r) where x is the axial distance from the 
input end of the tube and r represents the radial distance from the 
axis of tube symmetry. The four governing relations are the Navier- 


Stokes equations, the continuity equation, the equation of state for 
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an ideal gas, and the unsteady form of the energy equation. Under 
the assumptions that: 

a. The disturbances are small, 

b. The tube radius iS small compared to its length, and 

c. Laminar flow exists throughout the system, 
the controlling equations may be linearized in terms of the perturba- 


tions to yield the following set of coupled partial differential 


equations: 
AOD oot poeVie Malaga tt pose 
Bm + O- [at rel ee heat (A.02) 
op 
a ~ (A.03) 
ou , ov, V i: 
p 
H+ PH =p (A.05) 
O 
1 al : 
Nees 2 jwo = iwp ,gCoT (A.06) 


These partial differential equations are subject to the fol- 
lowing boundary conditions: 
a. At the wall of the tube, r = R; 
wix,R) = 0, vixgR) = 0,and TixsR) = 0. (A.07a) 
The latter boundary condition implies that the conductivity of 
the wall is assumed to be sufficiently large, therefore, the tempera- 
ture variation at the wall may be neglected. 
b. At the center of the tube, r = 0; 
by symmetry conditions, v(x,0) = 0, and 
u(x,0), T(x,0), p(x,0), and o(x,0) are finite. (A.07b) 
Equation (A.03) implies that the pressure perturbation is a 


function of axial distance x only. The energy equations may be treated 


20 


by a separation of variable technique, namely, let T = f(x)h(z) 


where: 
Z a By 
=—- ca 
R 
+ Pow, i 
« = 7? (FS) 


u 


Pr = Prandt] number = wis : 


Equation (A.06) becomes a form of Bessel's equation, namely: 








d?h(z) , 1 dh(z) i 
dz? " a dou h(z) = OogCp _— F(x) (A.08) 


which has a solution of the form: 


= i oe 
h(z) C.J (z) CoY (2) + 55gCP F(x) - (A.09) 


From equation (A.07b) it follows that C= 0. Sitee ater: 


T = 0, the function f(x) may be determined: 


: ! _ pix) 
f(x) aS C0 CaP ry 0 ,9CP q (A.10) 
Combining equations (A.09) and (A.10) yields an expression for the 
variation of the temperature perturbation amplitude as: 


Jot xP’) 


: i R p(x) 
T = f(x)h(z) = [1 - J,(ePr4) DogCP : (he) 
It should be noted that the Bessel] functions of equation 
3 : 
(A.11) have complex arguments due to the presence of i* =e iD 
the « constant; hence we introduce the Kelvin functions: 
3: 

Joly i7) = ber y + i bei y 

or alternatively in terms of amplitude and phase: ( Ags) 


2 


Joty 17) = Mo(y) [cose (y) + i sine (y)] . 


The first of the linearized Navier-Stokes equations, equation 


(A.02), can be similarly solved for the amplitude of the axial velocity 


5] 


perturbation, resulting in: 





J(5 a) 
Bee oe 1_ dp 
u = J(a) - | Twp, qx ° E181) 


Substitution of the temperature variation, equation (A.11), into the 


equation of state yields: 


0 = oo fi -¥ ri) - 5 “0 ae eda (A.14) 


Finally, the equation of continuity, equation (A.04), can be con- 
Sidered using the results for the axial velocity and density perturba- 
tions, equations (A.13) and (A.14), respectively; this becomes: 


v ips 
Ley) ae a , eye JOR nee. 
r or 1wo a ‘ Y J, aPr%) 


O 

J (ea) 

$2 lon f (A.15) 
dx? J Ca) 


Equation (A.15) may be integrated with respect to r in order to find 





(rev) with an arbitrary functional F(x) upon recognition that: 


atyd,(y)] = yoo(y) . 





Hence: 
J,(= aPr@) 
se ] Ww? 1 2 Y-] 1 2 rR R 
vt) m{iew ; baa 5 7 OPré Ta Pr) 
J, (Ga) 
d@P YR ~1‘R iL 2 
dx2 exer "6G "| a rf (A. 16) 


Applying the boundary condition, y(x,R) = 0, one obtains: 


. u)* R2 No Ae) Pr) R2 J ( ) d2 
nein A ee 2 Tay &e (A.17) 


uz 


where the Bessel function of the first kind and second order is 


obtained from the recurrence relation: 


Consideration of equation (A.16) in the limit when the radial coor- 
dinate, r, tends to zero requires that the constant F(x) equal zero, 
yielding the governing ordinary second order differential equation for 


the pressure perturbation: 


Jalen ) J Ale: ) 42 
az, YP I wh +s oe: 3 Tay ae * (Ae8) 


The .. equation provides a solution for the pressure perturba- 


tion: 
= A exp(ox) + B exp(-ox) (A.19) 
where 
ky N71 
,- Ps es ee veh y-1 92 loPré) 4 
ay J, (a) Y J (aPr’) 


The constants A and B may be determined after the boundary 
conditions at both ends of the tube have been prescribed, in addition 
to assumptions dealing with the instrument volume, namely: 

a. Pressure and density in the instrument volumes are 


only time dependent, and 


b. The pressure expansion in the instrument volume is a 


polytropic process, described by: 


P 





u c= constant . (A.20) 
(p,,) 


DS 


Consider the system as sketched below in a sectional view: 


Tube of length L 
radius R 








Instrument volume, 
use subscript "v" 





entrance to tube, exit of tube, 

use subscript "1" use subscript "2" 

If the length of the tube is L, then the pressure at the inlet, 
station 1, and at the exit, station 2, is given by equation (A.19) as: 


P, = A+B 


] 


Py = A exp(oL) + B exp(-oL) . (A.21a) 


The constants A and B can be expressed in terms of P, and P., as: 
-] 
A = a-sinnta_) LPy exp(-ol) - Pol 


B= serncacy Oy exp(oL) - P., ee (A.21b) 
The instrument volume has an influence upon the pressure response and 
is brought into consideration by requiring that the mass flow at 
Station 2, the tube exit, be in balance with the change in mass in 
the instrument cavity. The axial velocity amplitude at station 2 may 
be obtained by combining the pressure equation, equation (A.19), with 


the velocity solution, equation (A.13), to yield: 


J. (qe) 
u, = i, se - | [ exp(oL) - B ext -at | (A.22) 


At x = L, the mass flow can be expressed in a harmonic form by a 


straightforward integration process: 


54 


lwt lwt 
me.” Sgn (o,us(@nrdar 


meet = $2 R? tar exp(oL) - B exp(-oL) Jet : (A. 23) 
On the assumption that the pressure perturbation in the instrument 
cavity is equal to that at the tube exit, P> = Py» one obtains from 
a linear expansion of the polytropic relation of equation (A.20): 


Twt 





p Pe pee p 
ve © age V 7 6 
(F yk = constant “To +p elutyk soak 
Py O V Po 
which simplifies to: 
Pp = a, ~ Py (A.24) 


The mass in the instrument volume can be expressed in terms of the 
density under the assumption of constant volume. This assumption 
implies that the harmonic motion of the transducer diaphragm, which 


forms part of the cavity wall, is small. Hence: 





a iwt\ _ Y iwt 
Ly - v (eo, t fhe ) = V (eo, + az ok Poe ae (A.25) 
By conservation of mass, it is clear that: 
dm : 
= eee (A.26) 
Combining equations (A.26), (A.23), and (A.21), one obtains: 
J, (a) 
1WY a TO p2 2 ] 7 
a2 k “v P, * iw : J.(a) sinhCGLy [-Pi + p2 cosh(4L) ] 
or finally, the result for the dynamic response is obtained: 
Pe we ¥ Wh ts a) ; 
— = +=—T , 
7 [cosh(¢L) az ( vs 3 Tey Cue) (hee?) 


The formulation of equation (A.27) has been made by Bergh who 


has obtained good experimental confirmation. The frequency response is 


Significantly influenced by system geometry and mean pressure level, 
while the influence of the polytropic process index, k, was found by 
Bergh to be minor. In actual data reduction, one would take into 
account the value of mean static pressure in order to obtain a more 
precise definition of the dynamic response curve. This effect would 
be important in an experiment where the static pressure level might 
vary by more than about plus or minus 0.25 atmosphere. 

It is not difficult, as shown by Bergh, to extend these 
results to a series hookup of tubes and volumes by a repetitive 
analysis where the frequency response for input and exit of each 
individual tube length is obtained separately. It would be interesting 
to determine if a series connection of several tubes and volumes could 
lead to a more flat frequency response curve. This remains as an area 


for future investigation. 
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APPENDIX B 


PREDICTION OF UNSTEADY PRESSURES DUE TO 
HARMONICALLY OSCILLATING PLAIN FLAP 


The analysis described here may be considered as classical in 
the area of incompressible inviscid unsteady aerodynamics and is at- 
tributable to Kussner and Schwarz [6]. An introduction to the method 
is described in the aeroelasticity test by Fung [7], and for that 
reason only the particular terms involved in developing a solution 
for the pressures in the presence of an oscillating plain flap will 
be treated here. 

Consider a two-dimensional airfoil section of unit valued semi 
chord. A trigonometric coordinate transformation (sketched below) for 


chordal location may be used as is familiar in thin airfoil theory: 


Plain Flap 





> ».4 


Hinge 


Let x = -cosgf so that x = -] at the leading edge, 8 = 0, and x = +] at 
the trailing edge, 8 = 7. The hinge line of the plain flap is located 
by the angle wp. 

The harmonic motion of the airfoil is defined by: 

ies Oasis BiG) 
and 

Z = 6(cosp - cosp)e for p< B<m7. 


The downwash on the mean line can be described by: 


sy; 


or in trigonometric notation for the case of the oscillating plain 
flap: 

v(g,t) = use Lik( cosy - cosg) + 1] (B.01) 
where again py < 8 < m7 and the non-dimensional frequency term is: 


a Ge 
ko (ame Tale 


With this notation, a solution for the unsteady pressure developed at 
a chordal. location corresponding to the ‘angle 8 can be developed after 
considerable analysis. The results of Kussner and Schwarz can be sum- 
marized for the dimensionless pressure differential as: 

aCp(p,t) = + 6e'*[Re(k.6,y) + ilm(k,e,¥)1 (B.02) 
where 

Re (k,8,W) = Real part of expression 

Im (k,8,W) = Imaginary part of expression 


and 


Re (k.8.0) = CF(k) 0, (v) - © G(k) ®,(v) - simp] cot § 
+ L(8,¥)[1 + 2(cosy - cosg)] 
- KeL(n-)2 cosw+ sinw]sing -(r-W)sing cose 


+ L(8,W)(cosw - cosB)7} (B.03a) 
Im (k,8,¥) = [G(k)®,(v) + $ F(k)o,(v) - 55(W)Jeot $ 
+ 2k(1-W)sing. (B.03b) 
In equations (B.02a) and (B.02b), we further identify that: 
F(k) + G(k) = C(k) = Theodorsen coefficient which is tabulated 


in Reference [7] and: 


08 


_ ] 1-cos (8+) 
L(8.W) = 5 In Cae 
O,(y) = (m-v) + siny 
®,(p) = (m-y)(1+2 cosp) + sin (2+cosy) 


O,(y) = (m-y) + sin cosy 
A set of calculations were made to correspond to the wind-tunnel 
experiment, namely, the pressures at the 75 per cent chord station in 
the presence of an oscillating 20 per cent chord radius nosed plain 
flap. Results of the calculations relative to amplitude change and 
phase angle for the pressure differential are shown on Figure 16 where 
the values are normalized relative to the k = 0 (stationary flap ona 


fixed airfoil) situation. 
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APPENDIX C 


X-Z COORDINATES OF THE NACA 64,~-012 MODIFIED 
AIRFOIL USED IN THE WIND TUNNEL TESTS OF 
THE UNSTEADY PRESSURE MEASURING SYSTEM 
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